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ABSTRACT 



r- — Aims. The Two Micron All Sky Survey, along with the Stellar Population Synthesis Model of the Galaxy, developed in Besan^on, is used to 
^vj _ calculate the extinction distribution along diiferent lines of sight. By combining many lines of sight, the large scale distribution of interstellar 
^^ ' material can be deduced. 

^^ Methods. The Galaxy model is used to provide the intrinsic colour of stars and their probable distances, so that the near infrared colour 
^^ excess, and hence the extinction, may be calculated and its distance evaluated. Such a technique is dependent on the model used, however we 
>D I are able to show that moderate changes in the model parameters result in insignificant changes in the predicted extinction. 

Results. This technique has now been applied to over 64000 lines of sight, each separated by 15', in the inner Galaxy ( |/| <100°, \b\ <10° 
). We have projected our three dimensional results onto a two dimensional plane in order to compare them with existing two dimensional 
^ l' extinction maps and CO surveys. We find that although differences exist due to the different methods used or the medium traced, the same 
large scale structures are visible in each of the different maps. Using our extinction map, we have derived the main characteristics of the large 
scale structure of the dust distribution. The scale height of the interstellar matter is found to be 125^,^ pc . The dust distribution is found to be 
asymmetrically warped, in agreement with CO and HI observations of the ISM. However, the slope of the dust warp and the galactocentric 
^ , distance where it starts are found to be smaller than the values measured for the external HI disc: for positive longitudes the angle is 6 =89° 
;;• ■ , it starts at 8.7 kpc from the Galactic center and grows with a slope of 0.14 , while at negative longitudes, the angle of the maximum is at 
6 =272° , the starting radius 7.8 kpc and the slope 0. 11 . Finally, the presence of dust is detected in the Galactic bulge. It forms an elongated 
structure approximately 5.2 kpc long and lies at an angle of 30 ± 5° with respect to the Sun-Galactic centre direction. This may be inter- 
' $H preted as a dust lane along the Galactic bar. This resulting extinction map will be useful for studies of the inner Galaxy and its stellar populations. 
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1. Introduction Galaxy due to the amount of extinction and the high source 

confusion. The determination of interstellar extinction in this 



zone, which covers about 10% of the infrared sky, is thus an 
invaluable step to identifying low extinction windows at low 
Galactic latitude. 



Interstellar extinction remains a serious obstacle for the obser- 
vation of stars in the Milky Way, and for interpreting these ob- 
servations in terms of Galactic structure. It has two main ef- 
fects, that of dimming sources and reddening them. Any stel- 

lar study near the Galactic plane would need to coiTect for gg^^^^j ^^^^-^^ ^^^^ ^^^^ undertaken to further our knowl- 

this effect to enable a proper interpretation of the observations. ^^^^ ^f ^j^^ distribution of the extinction in our Galaxy. Most 

Recent infrared surveys provide a superb tool to probe further ^f ^^^^^ ^^^^-^^ ^ ^^^ dimensional map of the extinction in 

into the densest parts of the Galaxy as the extinction suffered in ^^^^^^^ ^^^-^^^ ^f ^j^^ q^I^^^ -^ ^j^^^j^ ^1^^^ calculate either 

theK, band is approximately one tenth of that suffered in the ^j^^ ^^^^ extinction or the total extinction for various hnes 

visiDie. q£ gjgjjf . Proebrich et al.l ( 120051) use accumulative star counts 

Furthermore, extragalactic studies often shun the so called ^^^^ ^j^^ 2MASS survey to construct a relative extinction map 

"Zone of Avoidance" at a round \b\ $10° and towards the mner ^f ^^e Galactic plane; .Sumj. C2004.) has created a two dimen- 

Send offprint requests to: D.J. Marshall sional map of the mean extinction towards the Galactic bulge 

* The three dimensional results will be available in electronic form using red clump giants from the Optical Gravitational Lensing 

at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79. 128.5) or Experiment II; .Dutra et aL C2003.) and .Schultheis et aL (.1999.) 

via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ determined the mean extinction towards the same region us- 
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ing 2M ASS data and DENIS data, respectively: ISchlegel et"ai] 
(119981) used the dust IR emission as observed by DIRBE to 
create an all sky two dimensional map of the total extinc- 
ti on integrated alon g eac h line of sight . Howe ver, according 
to lArce et all ( Il999l) . the ISchlegel et aP (Il998l) map overesti- 
mates the reddening by a factor of 1.3-1.5 in regions of smooth 
extinction with Ay > 0.5, and underestirnates it in regions of 
steep extinction gradients. ICambresv et alJ (120051) also noticed 
that the correlation of the near-infrared extinction to the far- 
infrared optical depth shows a discrepancy of about 30% for 
Ay >l mag, probably due to the ubiquitous presence of fluffy 
grains. 

Other authors have determined the distribution of extinc- 
tion along different lines of sight using various methods: 
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iFitzg erald' ('1968') provided the local distribution of extinction 
up to a few kiloparsecs using photoelectric photometry. Later 
on. lNeckel et alJ (Il980i) measured extinction and distances to- 
wards individual 0,B and A stars. These studies provide lim- 
ited resolution for the extinction distribution due to the scarce- 
ness of the distribution of hot stars, often related to clouds. 
They do not provide the distribution of the diffuse extinction. 
The study of Berdnikov et al. ( 1990) is limited to a sample of 
615 stars in 9 Kapteyn Selected Areas and shows considerable 
dispersion in extinction among the stars in a given direction. 
[Arenou et al. ( 1992) mapped the local extinction from a col- 
lection of spectral and photometric data. The validity of this 
map is limited to a distance of about 1 kpc, hence it cannot 
be used at \b\ <10°, where a significant part of the extinction 
is situated at distances larger than 1 kpc. In additio n, the reso- 
lution of the map is of the order of several degrees. ICambresvl 
( 1999a) determined detailed maps of extinction in nearby Giant 
Molecular Clouds from USNO star counts, which are very use- 
ful for studies of these clouds, but not outside. He also consid- 
ered the whole sky extinction from USNO star counts, but still 
with a two dimensional approach, and with a limited spatial 
resolution of 6° (C ambresv.l999 b). A great deal of other deter- 
minations can also be found in the literature which concern one 
or a few lines of sight. Putting all these maps together leads to 
inhomogeneous maps with poorly defined accuracy. They also 
have the problem of matching differen t results at the boundaries 
between authors dHakkila et alJl997l) . 

The relationship between infrared colour excess and ex- 
tincti on is well docurnented, and has b een in use for some 
time dJones et al."l984": !Lada et al.ll994 . In this type of study, 
stars that are reddened due to their location behind a cloud are 
compared to nearby, supposedly un-reddened, stars. Our ap- 
proach is to compare observed infrared colours with those from 
the Besangon model of the Galaxy. Using this technique, with 
which we are able to calculate the extinction along the line of 
sight, we have created a map of K, band extinction (Aks) in 
the plane of the inner Galaxy (|/| <100°, \b\ <15°) and in three 
dimensions. Our method for determining extinction does not 
assume any model for the dust distribution. Instead, we com- 
pare the difference in the J-Ks colour distr ibution of s tars from 
the Besangon model of the Galax y (Robin et al.l2003l) to those 
from 2MASS observations dCutri et al. 2003); the difference 
between the two is then attributed to interstellar extinction. The 
J-Ks colour index provides an excellent probe of interstellar 



extinction, in the Galactic plane where it is high, as radiation at 
two microns is less absorbed than visible light by around an or- 
der of magnitude. This technique can therefore provide a global 
picture of the three dimensional distribution of interstellar dust. 
Our use of the 2MASS data is discussed in ^ the gener- 
ation of catalogues of stars using the Besangon model of the 
Galaxy in ^ In ^we describe how we transform a difference 
in colour between observed and modelled stars into a prediction 
of the three dimensional extinction along a line of sight. The 
three dimensional results are presented in ^ these are com- 
pared to observations of CO and to other extinction data, and 
the large scale characteristics of its distribution are discussed. 
The sensitivity of the results to changes in model parameters is 
estimated and limits and possible bias are discussed in ^ Our 
results are summarised in ^ 

2. 2MASS Point Source Catalogue 

2.1. Description of the survey 

The Two Micron All Sky Survey (2MASS) is a ground based 
survey which uniformly scanned the entire sky in three near- 
infrared bands (7, H & Ks). Amongst its final products is the 
point source catalogue (PSC) which includes point sources 
brighter than about 1 mJy in each band, with signal-to-noise 
ratio (SNR) greater than 10, using a pixel size of 2.0". Each 
star has accurate pho tometry and astro metry as defined in the 
level 1 requirements dCutri et a l. 2003). 

The quality of the data was destined to meet or surpass all 
of the level 1 requirements. We list three of them below : 

- Photometric sensitivity of 10-cr at 15.8, 15.1, 14.3 mag at 
J, H, Ks respectively for \b\ >10° 

- Completeness >0.99 at 10-cr sensitivity limits 

- ReHabiHty >0.9995 

These requirements were met for the entire unconfused sky. 
However, the 2MASS survey has no level 1 requirements for 
the Galactic Plane - these specifications apply only to high 
Galactic latitudes. Within approximately a 6° radius of the 
Galactic centre, the average magnitude for a SNR=10 source 
is up to ~2 mag brighter than the typical SNR=10 source in 
unconfused regions. This is also seen for various regions in the 
Galactic plane at \b\ % 3° and + 60° from the Galactic centre. 

As well as containing default magnitudes and associated 
uncertainty for each star, the 2MASS PSC also contains a num- 
ber of flags to aid the user in selecting stars appropriate for 
the study. We make use of the Read Flag (rd_flg), which in- 
dicates the method used to determine the default magnitude 
of the source, and the Photometry Quality Flag (ph_flg) which 
gives an indication of the quality of the photometric determina- 
tion (signal to noise ratio, measurement error). More informa- 
tion on the differe nt flags and their meanings can be found in 
ICutrietal.l(l2003l) . 

The final 2MASS project also includes an extended source 
catalogue (XSC). By using a time var ying PSF as well as the 
J-Ks and H-K^ colour of the sources. lCutri et alJ (120031) were 
able to separate the extended sources from the point sources. 
After this classification, the final 2MASS PSC contains only 
point sources with a reliability of over 98%. 
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2.2. Completeness of the 2MASS PSC 

We are interested in regions in the Galactic Plane, where the 
completeness varies as a function of stellar density and can be 
as much as two magnitudes brighter than high latitude fields. 
The PSC can be compared quantitatively with model simula- 
tions only where the former is shown to be complete. As such, 
we must compute the completeness limit field by field. 

The 2MASS definition of completeness is "the faintest 
magnitude bin which recovers 99% of th e expected source 
counts". However, for their atlas images, ICutri et alJ ( l2003h 
suggest another method, described below, which we have 
adopted. Star count histograms are constructed for each field 
at 0.2 magnitude intervals for point sources with magnitudes 
between 9 and 18 mag and where a reliable estimate of the 
photometric error could be determined in the appropriate band. 
This latter condition excludes only a few percent of the total 
number of stars in a field except in high confusion fields, for 
example near the bulge or in the Galactic plane, where as many 
as 20% of the stars in the 2MASS PSC may not have an associ- 
ated photometric error. The completeness in a particular field is 
defined to be the bin that contains the highest star count. Stars 
detected by 2MASS but with photometry so poor that the error 
could not be determined are not used for the calculation of the 
limiting magnitude. This ensures that we do not overestimate 
the completeness of the observed stars which would introduce 
a large number of modelled stars with no observational coun- 
terpart. 



2.3. 2MASS selection criteria 

As stated, there are a number of flags for each star in the PSC 
in order to select them based on the quality of their photome- 
try, or to reject them due to the risk that they are contaminated 
by nearby bright objects. However, if we are too severe with 
our selection we will reject valid stellar detections. We will be 
comparing the selected stars with the Besangon model; it thus 
follows that we must select all point sources from the PSC that 
correspond to actual stellar detections. 

As we are interested in using the difference in observed 
and modelled J-Ks colour, we require detections in these two 
bands in order to admit a star into our calculation. However, 
a point source that has been detected in the J and Ks bands 
without a corresponding detection in the H band is not likely 
a valid stellar detection. Therefore, we restrict our selection to 
stars that have been detected in all three bands (2MASS rd_flag 
^ 0). Stars above the completeness limit in either the J, H or 
Ks bands are also rejected. We use the J - Ks colour index as 
it is less sensitive to metallicity and gravity than the H - Kg or 
J - H colour indices. 

This "full" catalogue includes all detected stars in the par- 
ticular 2MASS field, even those with less accurate photometry. 
However, as they are real detections they must be included in 
our comparison with the Besanfon model. 



2.4. Missing fields in the PSC 

The 2MASS catalogue contains a small number of "physical" 
or "effective" gaps, due either to unscanned areas of the sky or 
the removal of valid sources near tile boundaries, respectively. 
The total area of these gaps in our results amounts to less than 
0.04 square degrees. 

Other gaps in the 2MASS catalogue, which add up to a 
larger area in our results, include regions where bright stars 
have been removed, effectively creating a hole in the stellar 
density. All stars within the radius of the removed star are then 
excluded from the field of view; as these "missing" stars sam- 
ple all stellar populations and distances, they do not introduce 
any bias or error into our extinction determination. 

3. The Galactic model 

3.1. Description of the model 

The stellar populatio n synthesis model of the Galaxy con- 

J~~^ — II 1 

structed in Besan9on JRobin et alJI2003l) . hereafter called the 

Galactic model, is able to simulate the stellar content of the 
Galaxy by modelling four distinct stellar populations: the thin 
disc, the thick disc, the outer bulge and the spheroid. It also 
takes into account the dark halo and a diffuse component of 
the interstellar medium. It can be used to generate stellar cata- 
logues for any given direction, and returns information on each 
star such as magnitude, colour, and distance as well as kine- 
matics and other stellar parameters. 

The approach of the Galactic model is semi-empirical as it 
is based on theoretical grounds (for example stellar evolution, 
galactic evolution and galactic dynamics) and is constrained by 
empirical observations (the local luminosity function for ex- 
ample). The Galactic potential is calculated in order to self- 
consistently constrain the disc scale height. In addition, the 
model includes observational errors and Poisson noise to make 
it ideal for direct comparison with observations. 

The Galactic model has been developed to return results 
in the near-infrared and visible filters. It has been shown to 
reproduce the stellar content in the infrared to a high degree 
of accuracy (see ^6.2> . and is a powerful tool to extract the 
extinction information embedded in the 2MASS observations. 



3.2. Model parameters 

In this study, limited to the inner Galaxy, the bulge and thin 
disc are the dominant populations. Hence, the thick disc and 
spheroid model parameters will not be described here, but can 
be found in .Robin et aLL2003.) . 



3.2.1. Thin disc 

The thin disc is described by an evolutionary scheme with a 
two-slopes initial mass function and a constant star formation 
rate over the past 10 Gyr. The thin disc is divided into 7 age 
components, the first of which is called the young disc (age 
< 0.15 Gyr). The other six components define the old disc. 

The old thi n disc density distribution is modelled using the 
Einastd (119791) law. The distribution of each old disc compo- 
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nent is described by an axisymmetric ellipsoid with an axis ra- 
tio depending on the age; the density law is described by the 
subtraction of two ellipsoids : 



Pd - Pdo X [disc - hole] 



(1) 



where 

disc = exp 

hole — exp 
and : 



/0.25 + |A] J 2 



Rdl \exRd 



-.10.25 + I ^^ + 



RhJ \exRi 



- R and Z are the cylindrical galactocentric coordinates 

- e is the axis ratio of the ellipsoid. Recently revised axis 
ratios of the 6 age compon ents of the old thin disc can be 
found in'Robi netal] ( l2003l) : 

- 7?d is the scale length of the disc and is 2.510 kpc 

- /?h is the scale length of the hole and is 680 pc 

- the normalization p^o is deduced from the local luminosity 
function ( Jahreifi & Wielen 1997) , assuming that the Sun is 
located at R = 8.5 kpc and Z = 15 pc. 

The two values for R^ and Rh are recen t results from the 

I ^ " I ' ' 

method described in IPicaud & Robinf (2004). The equation of 

the old disc differs from that found in .Robin et al.. ^2003 ): the 

version presented here improves the modelling of the disc hole. 

The young thin disc density distribution is : 
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where a^ = 7?^ + (|) ,R,Z and e have the same significance as 
above and the scale lengths for the young disc are /?£/ = 5 kpc 
and Rh - 3 kpc. 

The stars from the young disc represent approximately 
0.5% of the stars in this study, thus the impact of this popu- 
lation in our method is negligible. 

3.2.2. Outer Galactic bulge 

An analysis of DENIS data in the region wit h Galactic coor- 
dinate s -8° < Z < 12° and \b\ < 4° allowed IPicaud & RobirJ 
( 1200 4) to constrain the shape of the triaxial bulge, as well as 
its age. The bulge is found to be boxy, prolate, with axis ratios 
1 : 0.30 : 0.25, and oriented 10.6 ± 3° with respect to the Sun- 
centre direction. A full description of the parameter values of 
the bulge density law can be found in Jicaud & RobiaC2004.) . 

3.2.3. Spiral structure and Galactic warp 

The version of the Galactic model we use here does not in- 
clude any spiral structure. The stars that we will be using for 
our method, mostly old K&M giants, are, according to the 
Galactic model, well over 1 Gyr old which corresponds to 
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Fig. 1. Grey-scale representation of the photometric error as a 
function of the ^5 magnitude for the 2MASS observations in 
the field (l,b) = (50°, 0°). The soHd line shows our best fit 
using Eq|3 



about 4 Galactic rotations. The effect of the spiral structure is 
thus assumed to be negligible for stars older than 1 Gyr. 

The stel lar warp is incorpor ated into the Galactic model as 
described in Robin et al. (^2003). Here we give just the relevant 
values. The Galactic warp is modelled as a shift of the galacto- 
centric coordinates perpendicular to the plane by a linear fac- 
tor equal to 0.18 beyond a galactocentric distance of 8.4 kpc 
(iDeiTiere and Robin 2001: Derriere 2001) . This effect is maxi- 
mum for an angle of 90° from the Sun Galactic centre direction. 

3.3. Diffuse extinction 

At the time of writing, the Galactic model available online 
(http://www.obs-besancon.fr/modele/modele.htinli 
is able to simulate interstellar extinction as diffuse extinction 
and as individual clouds. For the former, it is approximated 
using a double exponential disc with an ad hoc local normali- 
sation (usually 0.7 mag kpc ' for high latitude lines of sight). 
For the latter, any number of clouds can be added along the 
line of sight, defined by their distance and extinction (Ay). A 
cloud covers the entire field of view and therefore affects all 
the stars at distances greater than it. 

3.4. Pliotometric errors in tiie model 

The 2MASS PSC supplies photometric errors for each band 
and for each star, where the conditions (atmospheric, crowd- 
ing) permitted their determination. The Galactic model is able 
to simulate different photometric eiTor laws to emulate the 
response of different detection methods (CCD, photographic 
plates, etc). Photometric errors are included on each photomet- 
ric band independently, the error being a function of the mag- 
nitude (as usually observed). An exponential function is able to 
reproduce typical CCD photometric errors very well. 
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The error law, as implemented in the model, is : 

o-(m) = A + e*^^"-*' 



(3) 



where m is the magnitude in the band considered, and A, B and 
C are fitted to the set of data with which the model predictions 
are going to be compared. Fig[2shows the 2MASS photometric 
errors for the K, band in the field (Z, b)M5Qi°, 0°). The solid line 
is the best fit using Eq|3] 

The photometric errors in the 2MASS PSC are not constant 
across the entire sky, as they depend on local atmospheric con- 
ditions and stellar density. To perform the fit for a given field, 
we select the 2MASS stars in that field where a reliable es- 
timate of the photometric error could be determined in the ap- 
propriate band. The parameters of Eq.|3]are then determined by 
fitting the eiTor law to the selected 2MASS stars. This enables 
us to define an average error with respect to apparent magni- 
tude, field by field, for the stars in the 2MASS PSC. 

3.5. Limiting magnitudes in tiie model 

In order to closely model the 2MASS stars, we cut the mod- 
elled stars at the 2MASS completeness and at a J, H and Ks 
magnitude of 9 as a lower limit. The faint limit ranges from 
ihatj ~10.5 to ~14.4, and nij ~12.0 to ~15.8, depending on the 
crowding of the field. As the magnitude increases, the relia- 
bility of the 2MASS catalogue in the Galactic plane decreases. 
Therefore, to ensure the we have the highest reliability we fix 
the faint magnitude limit of the Ks band to a maximum of 12. 
The maximum magnitude used for the J band, however, is cho- 
sen to be the completeness hmit of the 2MASS observations in 
order to detect sources with high extinction. 

4. Extinction from 2I\/IASS and the Galaxy model 

Interstellar dust absorbs and scatters electromagnetic radiation, 
the effect being stronger for radiation at shorter wavelengths. 
This interstellar extinction results in observed stars appearing 
redder than they would in the absence of dust. 

The Galactic model provides us with simulated, initially 
unreddened stars and the 2MASS PSC contains observed, red- 
dened stars. If we assume that the Galactic model successfully 
predicts the average intrinsic colour and distance of stars in the 
Milky Way, then we may suppose that the difference between 
the modelled and observed stellar colour distributions is due 
solely to the effects of extinction. 

The method that we detail in this section enables us to de- 
termine the extinction, as a function of distance and along any 
particular line of sight, which minimises the difference between 
the modelled and observed stellar colour distributions. 

4.1. Colour distance relation 

In Fig|2we present an example field from the Galactic model 
towards 1=330°, b=1.5°; the J-K^ distribution is clearly a func- 
tion of distance. The dwarf and sub-giant populations are con- 
spicuous on the bottom left of Fig. |2 as they do not follow the 
same colour distance relation as the regular giant population. 
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Fig. 2. Simulated colour-distance distribution towards 
l=330°b=1.5°. Crosses are giants, triangles are dwarf stars and 
the diamonds denote rejected giants ( ^.3> . The group of stars 
in the lower left corner is dominated by dwarfs and sub-giants. 
There is a relation between colour and distance for the giants, 
first dominated by K&M stars, then, at larger distances by 
RGB/AGB stars. 

Dwarf stars are observed locally; they are too faint to be seen 
at large distances. As such they do not suffer much interstellar 
reddening, so the reddest (faintest) dwarfs detected will, con- 
trary to the giants, be the closest. These dwarf stars are removed 
from our stellar selection as described in ^4.3! Most of the sub- 
giants which are not on the same sequence as the giants are also 
removed in this process. 

The colour distance relation for the giants holds due to the 
fact that in the Galactic Plane the stars that dominate the counts 
at 9 < /:, <14 are mostly KO to K2 giants and RGB/AGB gi- 
ants; the latter dominate at larger distances due to their high 
luminosities. Stars which are homogeneous in intrinsic J-K, 
colour and in absolute magnitude Mk would have the same 
colour in a region with no extinction, showing a naiTOW his- 
togram of colour. As extinction is an increasing function of 
distance, the J-K, colour of these stars will naturally increase 
with distance. 

In the following we show how we compute the extinction, 
and its distribution in distance along the line of sight, for a par- 
ticular field. We use the Galactic model to simulate a variable 
proportion of K and M giants, as well as RGB/AGB stars, as 
a function of K magnitude along each line of sight, resulting 
in a realistic J-Ks distribution. In this way, the colour distance 
property of giants can be used to translate the shape of a J-K, 
histogram of 2MASS stars in a given direction into a distribu- 
tion of the extinction along the line of sight. 

4.2. NIR colour excess 

As mentioned above, interstellar dust has the effect of redden- 
ing starlight. Thus a difference in stellar colour (observed - 
intrinsic) can be used to deduce the amount of extinction the 
starlight has suffered. For a group of stars assumed to be at 
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a similar distance, the difference between their mean intrinsic 
colour and their mean observed colour (the mean colour ex- 
cess) yields the total mean extinction at that distance. Using a 
reddening law (for example Mathis 1990), we may write; 



Aks = 0.67 X E{J-Ks) 



(4) 



where Aj^s is the extinction in the band K^, E(J-Ks) is the 
J-Ks colour excess and the bar denotes that we have taken the 
mean for the group of stars. 

Instead of using the intrinsic colour of the stars as reference, 
we use simulated stellar colours from the Galactic model and 
we compute the difference in extinction between observed and 
modelled stars: 

6Ak, = 0.67 X [(7=^)obs - (7=^)siml (5) 

where the obs and sim subscripts denote observed and simu- 
lated colours, respectively. In other words, with Eq|3] we are 
calculating the extinction necessary to reduce, as much as pos- 
sible, the difference in colour between simulated and observed 
colour histograms. This equation may be used if the simulated 
stars have been reddened or not, as it simply translates the dif- 
ference in colour between simulated and observed stars into a 
difference in extinction. 

In order to calculate the extinction distribution as a func- 
tion of distance for a particular line of sight and field size, we 
proceed as follows. We start by applying a simple extinction 
distribution, as explained in ^3.31 to the simulated stars. The lo- 
cal normalisation is chosen in order to minimise the difference 
in the J-K^ colour distributions of the observed and simulated 
stars. The simulated, reddened stars are then cut at the faint 
magnitude limits for the field of the 2MASS observations. The 
dwarf stars, identified as explained below in ^4.31 are removed 
from both the observations and the model. 

Next, both the simulated and observed stars are sorted by 
increasing J-K, colour; due to the colour distance relation 
mentioned above, this means that we are effectively sorting the 
stars by increasing distance. We then bin the simulated stars 
by colour; the number of stars in a bin is chosen such that the 
median distance of the simulated stars in each successive bin 
increases. 

The observed stars are put into the same number of bins 
as the simulated stars, each bin containing the same relative 
number of stars: 

_ A^obs 
^obs, ~ ^ y '^mod, 



N, 



(6) 



mod 



where A^obs and A^mod are the total number of stars in the ob- 
servations and model, respectively, and «obs, and Mmod, denote 
the number of stars in a particular bin for the observations and 
the Galactic model, respectively. We may then obtain the mag- 
nitude of the extinction by applying Eq|5l bin by bin, and as- 
suming that the extinction calculated is at the median distance 
of the simulated stars in that bin. Using the relative number 
of stars makes the method less sensitive to a difference in the 
number of stars between the Galactic model and the observa- 
tions, hence to the assumed large scale structure of the stellar 
populations ( ^6.2> . 



ai 1 1 
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Fig. 3. Dwarf star removal for the field (/,/7)=(330,1.5). The 
triangles are dwarf stars, the diamonds are sub-giants and the 
crosses are regular giants. The lines running from top to bottom 
near J-Ks - 0.9 are the cut-off for dwarf stars before (dashed) 
and after(solid) the first determination of the extinction. The 
dotted line at niKs = 12 shows our faint magnitude cutoff. 



4.3. Removing the dwarf stars 

In Fig.|3lthe colour magnitude diagram of the Galactic model 
for the line of sight (/, b) = (330, 1.5) is shown. The giant and 
dwarf populations occupy distinct regions in the diagram: the 
dwarf population is on the left hand side, the giants are on the 
right. The line running down the middle indicates the upper 
limit of dwarf J-Ks colour as a function of magnitude (de- 
scribed below). We select all stars redder than this upper limit; 
what will be left will be almost exclusively K & M giants and 
RGB/AGB giants. 

For a given apparent magnitude, we are able to define a 
maximum J-K, colour for the dwarf stars as follows. If we 
suppose that, on average, the closest stars are at 100 pc, then 
we can assign this distance to the reddest dwarf stars as these 
will be the closest. The absolute magnitude of these dwarf stars 
at a given apparent magnitude can then be calculated, initially 
assuming that they do not suffer any extinction. Using the main 
sequence relation, we are able to convert this absolute magni- 
tude to a J-Ks colour index. This colour index thus defines 
the upper limit for the dwarf stars. However the dwarf stars do 
suffer some extinction; after an initial estimate of the extinc- 
tion has been obtained using the giants ( ^.2> . we redden the 
upper limit used according to the extinction found at the dis- 
tance of the dwarf stars. This gives us a new selection of giants 
with which we may recalculate the extinction. Fig|3lshows the 
J - Ks cut-off before (dashed line) and after (solid line) the first 
extinction determination. As can be seen, this is not a large shift 
as the dwarf stars are indeed local. 

4.4. Iterative method 

After determining the extinction distribution as described 
above we apply this extinction to the stars from the Galactic 
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model for the particular line of sight. To test the resulting dis- 
tribution we construct a J-K^ histogram of the adjusted model 
and the 2MASS observations (including the dwarf stars) - these 
can be compared by means of a;t'^ test. 

We comp ute the y^ statistic for the two binned sets of data 
as detailed in lPress et anill992h : 



X 



= Z 



( VA^obs/MnodWmod, - VA^mod/A^obs«obs,)^ 



(7) 



where «obs, (A^obs) and nmod, {Nmod) are the number of stars in 
the /* bin of the J-Ks histogram (total number of stars along 
the line of sight) of the observations and model, respectively. 
The term involving the total number of stars ensures that the;^'^ 
statistic will be lowest when the shapes of the two histograms 
are the same, regardless of any difference in the total number 
of stars. 

After having adjusted the extinction some stars will be be- 
yond the completeness limit for the field whereas others, in- 
trinsically bright, will be dim enough to be admitted into our 
selection. We therefore repeat the process using the adjusted 
stars and recompute the;^'^ statistic. We continue this iterative 
process until we find a minimum in the;\f^ statistic; the line of 
sight extinction giving rise to this value is then taken to be the 
best fit for the particular line of sight. This is usually accom- 
plished in under 10 iterations. 

This method supposes that the observed and simulated stars 
in a given bin are at similar distances. As we start with a simple 
distribution for the extinction this may not be the case, initially. 
However, each successive iteration improves the extinction es- 
timation and decreases the difference between the observed and 
modelled colour distributions. Assuming that the spatial den- 
sity of the giants is well modelled in the Galactic model, and 
that the extinction is uniform across the field ( ^6.3.2> . the as- 
sumption that the modelled and observed colour bins represent 
stars at similar distances is justifiable. 

As an example, the initial (after one iteration) and final re- 
sults for a sample field are displayed in Fig|3 The initial extinc- 
tion (diamonds) is significantly different from the final extinc- 
tion (asterisks), with nearly a 1 kpc difference in the location 
of the sharp rise in extinction at ~3 kpc. This difference results 
in a poorer fit to the observed colour distribution and therefore 
a large value for the X^ statistic (Eq0. 

We may summarise the various steps of the method as fol- 
lows: 

1 . Assume simple distribution for the extinction ^3.31 or apply 
extinction from last iteration 

2. Remove dwarf stars and cut observations and model at 
magnitude limits 

3. Sort both catalogues by increasing J-Ks colour 

4. Group simulated stars into colour bins such that the dis- 
tance of each bin > distance of previous bin 

5. Group observations into same number of bins as simulated 
star bins, each bin having the same relative number of stars 

6. Compute extinction correction using Eq|5l bin by bin 

7. Compute J-Ks histogram and;t'^ statistic 

8. Continue until the;^f^ statistic reaches a minimum 
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Fig. 4. Initial and final results of the extinction calculation for 
the field (I, b) - (345, 0). Top: Aks as a function of distance. 
Each diamond represents a bin from the initial calculation, the 
asterisks represent the bins in the final result and the error bars 
represent the mean absolute deviation of the extinction in the 
bin. Bottom: J-K^ colour distribution of 2MASS and adjusted 
Galactic model. The red solid line represents the 2MASS data, 
the bold solid line represents the final result and the broken line 
shows our initial calculation 



5. Results 

The full results are listed in Table 1, only avail- 
able in electronic form at the CDS via anony- 
mous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ Each row 
of Table 1 contains the information for one line of sight: 
Galactic coordinates along with the measured quantities for 
each bin {A^^, distance and respective uncertainties). All rows 
also contain the corresponding x^ statistic to indicate the 
quality of the fit between the modelled and observed J - Kg 
histograms. 

5.1. Results for individual lines of sight 

The results for three example lines of sight are presented in Fig. 
13 (/, b) = (57, 0), (330, 1 .5) and (5, 0). The top row presents the 
distribution of extinction along the line of sight; each asterisk 
represents one bin from our results and the error bar is the abso- 
lute mean deviation of the A^j from the median of all the mod- 
elled stars in the bin. Each simulated star at a given distance is 
attributed an extinction by performing a linear interpolation be- 
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Fig. 5. Results for the three fields at galactic coordinates (I, b) = (57, 0), (330, 1 .5) and (5, 0). Top: Aks vs. distance. Each asterisk 
represents a bin from our method, and the error bars are the mean absolute deviation from the median extinction in the bin. 
Bottom: Corresponding J-K, histogram of the model and the observations. The red line represents the 2MASS observations and 
the black line is the result of our method 



tween the two neighbouring bins. In order to test the extinction 
law thus determined, we compare the J-Ks histogram from the 
2MASS observations to the J-Ks histogram resulting from the 
application of our extinction to the Galactic model. The his- 
togram for each line of sight is presented in the bottom row of 
Fig-El As can be seen from these three examples, our method 
works very well in areas of high extinction. 



5.2. Two dimensional extinction maps 

By combining many lines of sight we are able to construct 
a three dimensional array of extinction values in the Galaxy. 
However, the values of A^j along each line of sight are not dis- 
tributed uniformly. If we want to know the extinction along a 
plane at a certain heliocentric distance we must first interpolate 
in between the values in order to get a regular grid. In what 
follows we have used linear interpolation between successive 
bins. 



5.2.1 . Comparison with two dimensional extinction 
maps 

Two dimensional exti nction maps towards th e Galactic cen- 
tre were c omputed by ISchulthei s et alJ (Il999l) (|/| < 8°, \b\ < 
1.5°) and iDutraetaTl 12003.) (\l\ < 5°, \b\ < 5°), hereafter 
S99 and D03 respectively. The method of S99 is based on 
the comparison of AGB isochrones to theoretical unreddened 
isochrones cBertelli et aL.1994.) . using DENIS data, whereas 
D03 compute the extinction by fitting the upper giant branch of 



(Ks, J-Ks) colour magnitude diagrams to a de-reddened upper 
giant branch mean locus of previously studied bulge fields. 

Although they used near-infrared data as we have, our ap- 
proach provides a different method for extracting the extinction 
information from the data. 

In order to compare our map with theirs, we first smooth 
their maps to our resolution of 15'x 15'. The S99 method yields 
results in Aks bu t they were p ublished in Ay', we use the ex- 
tinction law from lGlassl(ll999t) . Aks - 0.0S9Av, as they did to 
reconvert the Ay to A^s- The D03 map was published in Ak so 
we convert it to Aks using the conversion factor used in their 
article (A^ = Q.95Aks)- Finally, we integrate our extinction to 
the centre of the Galaxy, assumed to be at 8.5 kpc from the Sun. 

The results of this comparison are displayed in Fig. |6l 
which shows a comparison between our results and the S99 
map (left) and D03 map (right) at each pixel. 

Although the large scale structure of all three maps are very 
similar there is a systematic difference of ~ 0.2 - 0.3 A^s be- 
tween our map and that of S99. This is not surprising as they 
use the extinction law from G lassI lll999 ^. while we use the re- 
sults of Mathis. 1.1990.) . resulting in a different ratio for the K, 
band to J band extinction (A/cs/Aj). This ratio is used in the 
derivation of Eq|5land therefore has a direct impact on the ex- 
tinction determination. 

The extinction law in the infrared can differ significantly 
from one author to the next; the value we have adopted has 
been chosen so as to be consistent wi th the extinction law used 
for the Galactic model ( Mathisll990l). 

By combining the Mathis (1990) law, Ak/A i = 0.382 
with the relation Ak = 0.95Aks (.Dutra et al.ll2002l) we obtain 
Aks/Aj - 0.402, very close to the recent value (Aks/Aj - 
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Fig. 6. Co mparison between each pixel of our extinction map 
and that of'Schultheis e t al.l ( ll999l) (S99), for \l\ < 8°, \b\ < 1°, 
and Dutraetal. (2003) (D03), for \l\ < 5.0°, \b\ < 5.0°. Left; 
Aks of our map vs. that of S99. Right: Aks of our map vs. that 
of D03 . For both figures, the straight Hne represents the identity 
function. 



5.4. Three dimensional extinction 

In order to visualize the distribution of extinction elements in 
three dimensions, we divide the extinction between subsequent 
bins by the distance between them. The resulting map, in units 
of Ma: kpc"', gives us abetterideaof the location of obscuring 
dust along the line of sight. We assume a constant 6Ak kpc~' in 
between bins, as we have no information on the distribution of 
dust between bins. This results in an overestimation of the line 
of sight length of clouds for lines of sight with large distances 
between bins. 

These maps are presented in two forms: maps of the full 
area covered so far by our method at various distances (Fig.|8l 
and fully described in ^5.4. 1> . as well as a view of the Galactic 



plane from above (Fig.|9l referred to in ^5.4.21 



0.400) from llndebetouw et alJ (120051) . Howev er, this value 
is significantly different from loiassl (Il999l) (A/cs/Aj - 
0.347), which is close r to the value recently measured by 
iNishivama et al.1 JlOOdi . A^s/Aj - 0.331, in their study of the 
Galactic bulge. 

The application of the lMathis ^("l990") law to the S99 results 
effectively removes the systematic shift seen in Fig|6l D03 al- 
ready use the same value foT Aks/Aj as we do; the correlation 
between our map and theirs is very good. The slight depar- 
ture from the identity function at low extinction, comes from 
the poor sensitivity of the J-Kg colour index at these relatively 
low extinctions ( ^6.4. 1> . These low extinctions do not exist in 
the S99 map as the area covered by their map does not rise 
much above \b\ ~ 1.5; in addition to the more central regions 
of the Galaxy, the D03 map also probes higher latitude, lower 
extinction lines of sight. 

5.3. Comparison with velocity integrated CO emission 

IPame et alJ (120011) have created a large scale CO survey of 
the entire Milky Way, using the rotational transition 1 -0 of the 
CO molecule at 115 GHz, from new observations and existing 
CO surveys. By comparing the low velocity (local) CO emis- 
sion with an optical panorama of the inner Galaxy, they con- 
clude that there are few dark interstellar clouds where CO is 
not present or where the temperature is so cold that the CO 
molecule fails to radiate at 115 GHz. Galactic CO emission 
traces only dense regions of the ISM whereas our map also 
traces the diffuse component. Nevertheless we expect that both 
maps should trace the large scale distribution of interstellar 
matter in our Galaxy. 

Fig. shows the total extinction detected for each line of 
sight usi ng our method as w ell as the distribution of CO as pre- 
sented in lPame et alJ (1200 ih . where the two maps overlap. The 
resolution of the extinction map is 15'x 15'; the resolution of 
the CO map varies from 9' to 30', the region that we present 
here is mostly at the higher resolution of 9'. The main struc- 
tures that appear in both maps are indicated. The Musca Dark 
cloud, (/. b) = (-5 9, -9), that appears in our map is not in the 
IPame et alJ ( 120011) map as this region was not observed in their 
composite survey. 



5.4.1. Distance slices 

The three dimensional extinction towards the inner Galaxy is 
presented in Fig. |8] Each image shows the average local ex- 
tinction at distance intervals of 1 kpc. The scale is the same 
for each image and runs from to 0.3 A^-j kpc '. The solid 
line superimposed at each distance shows the stellar warp as 
implemented in the Galactic model. 

One obvious feature in the images is the geometrical effect 
of distance; high latitude extinction is mostly local so the disc 
is seen to shrink as we go further away from the Sun. In the 
first couple of slices, the local features mentioned earlier are 
visible: the Aquila Rift (30°, 4°) and the Local Arm towards 
Cygnus X (80°,0°). 

The giant molecular cloud (GMC), and supernova remnant, 
W44 is visible in the second and third slices (2kpc & 3kpc) at 
(35°, 0°). At distances between 3 and 5 kpc the tangent to the 
Carina arm can be seen at 78°. 

After 6 kpc, a lack of absorbing materi al is noticeable to- 
wards the Galactic bulge, an effect noted bv iGlass et al.l (119871) 
who found that the bulge has little dust compared to its CO 
emission. At these distances another feature of the disc which 
starts to appear is the Galactic warp. It becomes more pro- 
nounced at large absolute longitudes as we move farther out 
from the Sun. 

Finally, the lack of line of sight resolution for certain fields 
is obvious at (/, b) ~ (30, 5) and at distances greater than 3 kpc. 
Due to the fact that we use linear interpolation between suc- 
cessive bins and therefore do not assume any structure for the 
distribution of the dust in between bins, the large distance be- 
tween successive bins creates an unrealistically long structure 
along certain lines of sight. The extinction estimation itself is 
still assumed to be correct, within error. 

5.4.2. Tine Galactic plane 

By creating a map showing the local extinction in SAks kpc \ 
we are able to visualise the location of absorbing matter in the 
inner Galaxy. A view "from above" is displayed in Fig.|9lwhich 
is the average 6Aks for \b\ < 0.25 towards the inner Galaxy. The 
Sun is located at the middle left of the graph, and the Galactic 
centre is marked by an 'X'. 
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Fig. 7. Top: Total extinction integrated along the line of sight. Bottom: CO velocity integrated spatial map bv lPame et alJ ( l200ll) . 
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Fig. 9. Location of absorbing dust in the Galactic Plane (\b\ < 0.25). The Galactic centre, assumed to be at 8.5 kpc from the Sun, 
is marked with an 'X'; the Sun is located in the middle on the left side of the plot. The dotted lines are placed every 5 kpc from 
the Sun. Ticks on the border indicate Galactic longitude. 
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Fig. 10. Variation of the density of absorbing matter as a func- 
tion of distance from the Galactic plane. 



Many structures are visible in this figure, such as the lo- 
cal arm extending outwards from the Sun at / ~80°. Two tan- 
gents to the Sagittarius-Carina arm are also visible, extending 
from the Sun in two directions at I ~ 280° and I ~ 50°. The 
outer Perseus arm is above the top edge of the map and is 
therefore not readily identifiable. The inter-arm region between 
the Centaurus-Crux and Norma arms is visible at ~ 4.0 kpc 
from the Galactic Centre. After the Norma arm, centred on the 
Galactic centre is a large void of dust apart from an elongated 
structure which may well indicate the presence of dust in the 
Galactic bar, which we discuss below. 
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Fig. 11. Local extinction at 8 kpc from the Sun for positive lon- 
gitudes (top) and negative longitudes (bottom). The stellar warp 
as modelled in the Galactic model is shown as a dashed line; 
the solid line shows the best fit to our extinction data ( ^5.5.2> . 
The colour scale is the same as in Fig|8l 



In order to measure the scale length of the detected inter- 
stellar matter, we tested two profiles : 



p = posech (z/h^) 



(8) 



5.5. Galactic parameters 

We concentrate on determining large scale characteristics of the 
3D dust distribution from our 3D extinction map. The mean 
scale height, as well as the shape of the warp are measured. 
An estimate of the central dust bar orientation and size is also 
given. 



5.5.1 . ISM disc scale height 

Our 3D extinction map allows us to calculate the scale height 
of the dust component of the ISM The extinction between two 
successive bins was c onverted to a hydr ogen column density 
using the results from BohUn et alJ ( ll978h : this column density 
was then divided by the distance between the bins to obtain a 
hydrogen density per volume. This supposes t hat the metalhc- 
ity is t he same as in the solar neighbourhood: IChiappini et alJ 
( 1200 11) show that the abundance gradient in the Milky Way is 
relatively flat between 4 and 10 kpc from the Galactic centre 
but that it steepens sharply in the outer disc. To minimise the 
effect of any such gradient on our calculation, we avoid all ar- 
eas within 4 kpc of the Galactic centre and beyond the solar 
circle. We then calculated the mean hydrogen density at every 
10 pc from the Galactic plane. 



p = poe 



i-Wh:) 



(9) 



where /;, is the scale height of the ISM and po the hydrogen 
density in the Galactic plane. The use of the sech^ profile re- 
sults in a better fit and is shown in Fig^| The exponential 
gives /ij = 134;^|j pc andpo - 2.5 ±1.3 atoms cm"-', whereas 
the sech^ pr ofile yields hy = 1 25"^'^ pc and pn - 1.8+1 .0 atoms 
cm^l Both iMalhotral ll 19951) and Na kanishi & Sofuel ll2003h 
find that the scale height of the HI disc in the Milky Way 
varies as a function of distance; for both the mean value be- 
tween 4 an d 8 kpc is of the ord er of ~ 160 pc. The dust model 
of .Drimmel and Spergell ( 120011) . adjusted to the FIR emission, 
includes a scale height of 134.4 + 8.5 pc. We were not able 
to detect any significant flaring of the interstellar matter disc 
within the Galactocentric radii we tested. 

5.5.2. Galactic warp 

The Galactic warp is obvious in figure |8] especiaUy at large 
absolute longitudes and heliocentric distance. The magnitude 
of the warp in the Galactic model is calculated using Eg. 1101: 

Zwarp = r * (^ - '"warp) * C0S(6I - e^waip) (10) 

where Zwarp is the perpendicular shift from the mean ^ = 0° 
plane, R, and z are Galactocentric coordinates, y is the slope 
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of the displacement of the mid plane with respect to the plane 
of symmetry, r„aip is the radius after which the warp starts and 
^warp is the angle of maximum warp. Using our results, we are 
able to locate the z with the highest density for each (r, 9) by 
calculating the best fit with Eq[ro| We find that the warp is 
not symmetrical; for positive longitudes we find y - 0.14, 
=89° and rwaro = 8.7kpc and for negative longitudes 



-'warp 



' waip 



we find y - 0.11, 0warp =272° and rwaip = 7.8kpc. FigfTTI 
shows the local extinction per kpc at a heliocentric distance 
of 8 kpc, overlaid with the warp parameters of the Galactic 
model (dashed line) and with the above parameters (solid line). 
These values correspond to an absolute vertical shift of the 
plane al R - 12 kpc of 460 pc at both negative and positive 
longitudes, significantly smaller than the modelled stellar warp 
which predicts a vertical shift of ~ 650 pc at /? = 12 kpc. These 
new values for the dust warp are slightly higher than those of 
brimmel and Spergell (l200ll) (y = 0.0728, rwaip = 6.993 kpc), 
who predict a maximum vertical shift of 364 pc at /? = 12 kpc. 
This shi ft is different for the ste l lar wa rp, which they find to be 
higher. iLopez-Corredoira et al.l ( 120021) use a different formula- 
tion to model the warp and find the maximum vertical shift of 
the stellar warp aXR - 12 kp c to be 547 pc, in good agreement 
wit h the gas warp foun d by BurtonI (119881) . 

iNakanishi & Sofuel (.200 3 ) find that the Galactic HI disc is 
warped and asymmetric, with a larger shift at positive longi- 
tudes, as we do. They report that the warp is at a maximum for 
e = 80° and 6 = 260°, and that the warp starts at /? = 12 kpc, 
contrary to our values at around R - % kpc. At positive lon- 
gitudes the warp continues to reach a maximum displacement 
from the b - Q° plane of 1.5 kpc aXR - 16 kpc, resulting in a 
slope of 0.375. At negative longitudes they find that the warp 
is less severe, resulting in a shift of 1 .0 kpc below the b - Q° 
plane ?AR - 16 kpc, a slope of 0.25. Their values for the gradi- 
ent of the warp are significantly higher than ours, however our 
map does not have information out to these large galactocen- 
tric distances and so we are not able to constrain the large scale 
attributes of the int erstellar warp in the outer Galaxy. 

Most recently iLevine et al.l (12006 ) have shown that the 
outer HI warp can be described by a superposition of the three 
low mode (m = 0, 1 and 2) vertical harmonics of the disc, one of 
which (m = 1 corresponding to the integral-sign shape) dom- 
inates the warp for R 'i \5 kpc. It seems that it is this mode 
that we see in the dust from our study, which does not extend 
far enough into the external regions of the Galaxy to detect the 
other modes. They find that this mode is approximately linear 
for 10 :$ /? :S 15 kpc and that it has a slope of 0.197, signifi- 
cantly higher than our value. 

From these comparisons one sees that the warp is present in 
all galactic components (dust, gas and stars), all with the same 
node position and all asymmetric. However, the amplitude of 
the warp seems to depend slightly on the component one looks 
at: the dust warp (our study; Drimmel and Spergel 2001 ) seems 
to be less pronounced th at the ste llar warp (the Galac tic model; 
iLopez-Corredoira et alJ.2002: .Drimmel an d Spergel 2001), it- 
self less pronounced tha n the HI warp (Nakanishi & Sofua 
I2003r ILevine et al . 2006). However, this conclusion has to be 
taken with caution because the error bars are still large and our 
study does not go far enough into the external Galaxy. 




Fig. 12. Straight line fit that minimises the mean absolute devi- 
ation of the points of highest density along the Galactic bar. 



5.5.3. Dust bar 

The visualisation of our extinction map in the plane Figure- 
veals an elongated structure passing through the Galactic cen- 
tre. To determine its parameters, we fit the straight line to this 
bar that minimises the mean absolute deviation of the points of 
highest density along the bar (Figll2>. 

We find that the dust bar is at an angle of 4> =30 ±5° 
relative to the Sun - Galactic centre direction. The structure 
thus defined has a length of 5.2 kpc . The angle found is 
higher than the modelled stellar bar ( ^3.2.2> and that found 
by Babusiaux & Gilmore (2005) (0 = 22.5 + 5.5°), but still 
within the uncertainty of th e latter. Further more it is compati- 
ble with the values cited in iGerhardI ( 1200 ll) . who find that the 
angle may lie between ^ — 15° and (j) - 35°, although the bar 
angles based on stellar s tudies are systematically lo wer than 
the value we have found. iBissantz & Gerhard I ( l2002h find that 
COBE/DIRBE observations are compatible with an angle for 
the stellar bar between <^ = 15° and = 30°, but that the model 
with an angle of (p - 20° fits the observations best. These re- 
sults for the stellar bar suggest that the dust bar we have de- 
tected may in fact be dust lanes which precede the stellar bar 
at negative longitudes and trail it at positive longitudes, a fea- 
ture of our Galaxy suggested by Calbet et al. ( 1996). However, 
we detect a difference in extinction between positive and nega- 
tive longitudes of ~ 0.1 Aks, much smaller than their value of 
between 1 and 2.6 Ak- 

5.6. Comparison with other 3D methods 

As detailed in ^ other methods have been developed to mea- 
sure interstellar extinction along the line of sight. We now 
present our results alo i igside tho se of Neckel et al. (1980) , 
iLopez-Corredoira et alJ (l2002h and lOrimmel et al.l ( l2003l) . 

The three methods are described briefly below; for further 
information, please refer to the corresponding article. 

iNeckel et alJ (Il980l) : Extinction values and distances are given 
for over 1 1000 O-F stars. The intrinsic colour and absolute 
magnitude of the stars are deduced from the MK spectral 
type and from the Hfi index using the relationship between 
Mv and yS. 
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Fig. 13. The extinction, along the line of sight , in six different fields. For each field, we have compared our results (a sterisks) with 
the results of lNeckel et alJ (1 19801) (triangles'). IPrimmel et alJ ( 120031) (solid line) and lLopez-Corredoira et alJ ( 120021) (diamonds), 
where available. 



iLopez-Corredoira et alJ ( 120021) : This method is based on 
2MASS observations of red clump giants, for which a value 
of the absolute magnitude and J-Ks colour is assumed 
{Mk, = -1.65 and {J-K,\) = 0.75). (J-K,,MkJ CMDs 
of the giants are built for 1° x 1° fields and their displace- 
ment from an unreddened CMD, along with their absolute 
magnitude, are used to calculate the extinction and thus the 
distance to the extinction. 

Drimmel et al . (2003): The three dimensional dust model by 
IPrimmel and SpergeL ^2001.) is used to calculate the three 
dimensional extinction. The dust model is adjusted using 
the far-infrared residuals between the DIRBE 240 fiva data 
and the emission predicted from the model. 

Results for va rious lines of sight are presented in Fig. [O] 
The INeckel et al.l tl980.) data is the one with the largest dis- 
persion: young and hot OB stars are often surrounded by their 
parent cloud and are thus subject to possibly large and vary- 
ing circumstellar extinction and theref ore not representative of 
the wh ole field. Our results follow the iLopez-Corredoira et alJ 
(I2OO2I) results quite closely; this is to be expected as they are 
both based on the near infrared colour difference of giant stars. 
However significant dififerences do exist along certain lines of 
sight; this may be due to the fact that they do not use a model to 
determine the distance but instead use a unique absolute mag- 
nitude for all stars. Furthermore they do not use the same stel- 
lar populations as us, isolating the K giants from dwarfs and 
M giants using the "SKY" model dWainscoat et alJll992l) . In 
areas where this single population hypothesis doesn't hold dif- 
ferences between our results would be expected 

Fin ally, we are not always in agreement with lDrimmel et alJ 
(I2OO3I) . This is also to be expected given the region under study; 



anomalous dust temperatures, as well as lines of sight passing 
through spiral arm tangents or n ear the Galactic centre may in- 
troduce systematic errors in the IPrimmel et al.l (120031) extinc- 
tion map. Contrary to theirs, our map is not sensitive to dust 
temperature or to an assumed spiral structure but to the mod- 
elled K&M giant luminosity function (see ^6.2.2t . 



6. Discussion 

6.1. Use of the extinction maps 

Our results are available in the form of an electronic three di- 
mensional map. The map is presented in A^s, as the extinction 
law in this wavelength range is nearly universal ( Mat hisI 19901) 
and so does not depend on the ratio of total to selective extinc- 
tion. By adopting a s uitable extincti on l aw (fo r exam ple lMathisI 
il990), Rieke and L ebofskvl ( 1 19851) or lOlassI (Il999l) ). one can 
convert this K^ band extinction to an extinction in the visible 
portion of the spectrum or any other 

It was not possible to apply our method to some lines of 
sight in the region we have studied (|/| <100°, \b\ <10°) ei- 
ther because there were an insufficient number of stars in the 
2MASS observations for our pixel size or the extinction was 

In these cases we 



below the sensitivity of our method ( ^6.4.1 
have substituted the results of our method with the diffuse ex- 
tinction disc of the Galactic model ( ^3.3> . with a local normal- 
isation chosen to minimise the x^ test described in ^4.41 This 
solution was also preferred in regions where our method could 
not improve on the extinction predicted by the default extinc- 
tion of the Galactic model. To aid the user of our map, we have 
included details on each line of sight such as the;i'^ statistic and 
magnitude limits of the 2MASS catalogue. 
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Our results are most reliable in fields with homogeneous 
extinction. As this is not always the case, special care should 
be taken when using fields with large values for the;if^ statistic. 
In these cases, our method was not able to find a satisfactory 
solution for the line of sight extinction. 



6.2. Effect of changing the model parameters 

The determination of extinction using the Galactic model can 
be affected by changing the model parameters. The magnitude 
of the extinction that we calculate depends on the difference 
in colour between observed and modelled stars; the distance 
to the modelled stars (and hence to the extinction) is mostly 
dependent on the assumed absolute magnitude of the dominant 
population. 

This feature of our method means that the number of stars 
does not affect the extinction determination, so that if there are 
a number of stars in the model which are not present in the 
observations, we do not "invent" extinction in order to remove 
the excess stars from the model. 

The performance of the Galactic model to reproduce ob- 
servations in low extinction windows is a useful m easure of the 
validit y of the model. This comparison was done in Robi n et alJ 
(I2OO3I) : star count predictions (stars per magnitude and per 
square degree) in the Ks band were in good agreement with 
observations from the DENIS survey and so we refer the inter- 
ested reader to the aforementioned publication. 

We now briefly describe the effect on the extinction deter- 
mination of: the size of the hole in the modelled disc distribu- 
tion ( ^6.2. 1> : changing the bulge luminosity function ( ^6.2.2> ; 



and changing the stellar density ^6.2.3l on the resulting extinc- 
tion distribution. 



6.2.1 . Galactic disc hole 

In FigOwe show the median difference between each pixel 
of two extinction maps (mapl minus map2), at different helio- 
centric distances and for \l\ < 10° and \b\ < 2.5°. The Galactic 
model used for mapl has a central hole in the disc population 
of 1360 pc (R;,=1360 in Eq. [0 and the model used for map2 
includes a hole with a diameter of 680 pc. The scale in the 
figure runs from -0.1 to 0.1 Aks', a positive difference indicat- 
ing more extinction in the large hole model. As can be seen 
in Fig ll4l for distances closer than 6 kpc the large hole model 
predicts more extinction. We use the Galactic model to calcu- 
late the distance to the extinction: with fewer "distant" stars in 
the big hole model, we place the detected extinction closer than 
for the small hole model. The median difference is never large, 
however, and stays within the error of the method(' ^6.4.1> . 



6.2.2. Bulge luminosity function 

The effect of changing the bulge luminosity function on the de- 
termination of the extinction has been tested. Three luminosity 
functions were used : 
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Fig. 14. Influence of the assumed disc density law near the 
Galactic centre. Median extinction difference, at different he- 
liocentric distances and at \l\ < 10° and \b\ < 2.5°, between 
two determinations of the extinction using models with dif- 
ferent values for the disc hole scale length {R/, in Eq[2l. The 
error bars represent the mean absolute deviation from the me- 
dian. The difference is calculated per pixel, at each heliocentric 
distance. A positive difference shows that the model with the 
larger /?/, predicts a larger value for the extinction. 

Padova 7.9 & 10.0 Gyr Deduced from theo retical isochrones 

by the Padova team fcirardi et alJl2002l) . of which we test 

two bulge ages: 7.9 Gyr, 10 Gyr 
Bruzual 10.0 Gyr taken from the e volutionary bulge sy nthe- 

sis models of Bruzual & Chariot ( iBruzual et alJll997h . We 

test their 10 Gyr model. 

all of them based on a Salpeter initial mass function (a - 2.35) 
and assuming a single epoch of formation (starburst) as well as 
a mean solar metallicity (Z » 0.02). 

Although these luminosity functions alter the distribution 
of extinction along certain lines of sight, the overall effect is 
minimal. The extinction resulting from the use of each lumi- 
nosity function was compared at intervals of 0.5 kpc for var- 
ious lines of sight towards the outer Galactic bulge and the 
maximum difference in Aks between the resulting extinction 
distributions at each distance point was recorded. The ensuing 
difference distribution has a mean of 0.021 A/f^ and a standard 
deviation of 0.018 A/^^. As this value is comparable to the error 



of our method ( ^6.4. H . we choose to only use the default lumi- 
nosity function for the Galactic model, namely that of Bruzual 
10.0 Gyr 

6.2.3. Galactic scale length 

When modelling the Galactic disc, the density is normalised 
at the solar position. Therefore, in the plane of the Galaxy, the 
stellar density of K & M giants is largely determined by the 
scale length of the old disc (Eq.0. As the scale length shortens, 
the number of stars towards the Galactic centre increases, while 
an increase in the scale length has the opposite effect. 
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Fig. 15. Influence of the assumed disc scale length. Simulated observations are computed with a scale length of 2.51 kpc (first 
column), 2. kpc (2nd col.), and 3. kpc (3rd col.) while the scale length of the model used for computing the extinction is 
assumed to be 2.51 kpc in each case. Top: the distribution of extinction along the line of sight at (Z, fe) = (10°, 0°). The solid line 
represents the known extinction law of the simulated observations and the asterisks are the results of our method. In each case the 
scale length used in the simulated observations is included in the graph. Bottom: the coiTesponding J-K^ histograms with the 
extinction applied. It is noticeable that the fact that the model does not reproduce the number density of stars in the observations 
does not influence the determined extinction. 



In order to test the robustness of our method to a change 
in stellar density, we replace the 2MASS observations in our 
method with simulated observations with a known extinction 
distribution using the Galactic model. In the following, we vary 
the Galactic scale length in the simulated observations. We then 
proceed as normal (@, using the default model to determine 
the extinction. This simulates the impact on our method of 
choosing a value for the scale length of the Galaxy which is 
up to 25% different from the observations. 

In Fig[2]we show the effect of setting the scale length of 
the simulated observations to 2510 pc (the default value for 
the Galactic model), 2000 pc and 3000pc. The top row of this 
figure shows the extinction distribution of the simulated obser- 
vations (solid line) and the extinction distribution found using 
our model (asterisks with error bars). The simulated observa- 
tions have been subjected to a diff'use extinction as described 
in ^3.3l with a local normalisation of 0.7 mag kpc"' plus two 
clouds at 2 and 6 kpc. The bottom row shows the J-K^ his- 
togram of the simulated observations and of the default model 
with the calculated extinction distribution. 

The calculated extinction is found to reproduce the simu- 
lated observations' extinction very well, and always within the 
error bars. At large distances a small systematic difference is 
seen : if we assume a value of 2510 pc for the scale length 
while setting it to 2000 pc in the simulated observations, we 
produce an underestimation of the extinction of ^ 10% at dis- 
tances larger than 12 kpc. Underestimating the scale length (as- 
suming 2510 pc when it is 3000 pc in the observations) results 



in a slight overestimation of the extinction at large distances; at 
shorter distances, the position of the first cloud is also displaced 
by 0.5 kpc. 

6.3. Bias and uncertainty 

6.3.1. Circumstellar dust 

As we are using K&M giants as a tracer of interstellar ex- 
tinction, our stellar selection will be contaminated by AGB 
stars with circumstellar dust, making them redder than can be 
accounted for by interstellar extinction alone. M giants with 
spectral type later than M5 also show signs of circumstellar 
dust, although the amount has been shown to be negligible 
( iGlass and Schultheis..2003.) . In most parts of the Galaxy the 
ratio of AGB stars to RGB stars is small, of the order of a few 
%. Only in the inner parts of the Galactic bulge (at |/| < 2 and 
\b\ < 2) does the fraction of AGB stars become significant (up 
to ~ 20%). However, number densities of the different types of 
AGB stars such as semi-regular variables, Mira-type variables 
and the high Iti minosity OH/IR stars are still very poorly known 
(lH ahingll99fil). 

Ijura and KleinmannI ( 119891) studied the number densi- 
ties of mass-losing AGB stars in the solar neighbourhood 
and estimated t heir su rface density at around ~ 25 kpc"^. 
iLe Bertre et alJ (120031) used the Japanese space experiment 
IRTS to investigate the spatial distribution of mass-losing AGB 
stars. They found that objects with high mass-loss in the or- 
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der of 10"^ < M^lyr < \Q^^ dominate the replenishment of 
the ISM, however these sources constitute a minority of the to- 
tal number of stars in our sample (~ 10%). These stars have a 
large near-IR excess in J-K (> 0.5 mag). 

We conclude that there is only a very small level of contam- 
ination by high-mass losing AGB stars in our selection, mainly 
in the region of the inner bulge. Our values for the extinction in 
this region may be slightly overestimated at large distances. 

6.3.2. Cloud substructure 

For each line of sight measurement, the extinction that we de- 
tect is ass umed to cover the entire field at the given resolution. 
However, iFaison and Goss I (1200 ll) have resolved structure in 
the neutral interstellar medium using the Very Long Baseline 
Array at scales as small as 10 AU. The resolution of our map is 
15'; as we cannot resolve structure closer than ~ Ikpc, we will 
not resolve any cloud smaller than ~5 pc across. 

Nevertheless, as we are placing ourselves in the context 
of the large scale fluctuations in the ISM, the effects of these 
small clouds do not have a significant effect on our estimation. 
Clouds with an angular size much smaller than our pixel size 
will introduce anomalous reddening for a small fraction of the 
stars in a given field. Our use of the median as an estimator 
excludes these outlier points. 

6.4. Limits 

6.4.1. Maximum and minimum extinction detectable 

At low Aks, stars that are detected in the K^ band are almost cer- 
tainly detected in the J band. However, as column density in- 
creases and the reddening of the sources gets larger the chances 
of detecting a source in the Ks band without a corresponding 
detection in the J band increases. 

No selection effect is introduced, however, as we reject the 
modelled stars in the same way as they are in the data (same 
limiting magnitude in the J and K, bands). Instead, this will 
introduce a maximum Aks in the map as very red stars will not 
be detected. 

This upper limit depends on the magnitude limit in the / 
band of the 2MASS observations. As we place a lower bound 
on the Ks magnitude of 9, the highest J-K^ colour index ob- 
servable along a particular line of sight would be 7corap - ^.s(min) 
where the subscript comp signifies the completeness limit and 
min the brightest magnitude used due to our self imposed limit. 
The highest extinction observable along a line of sight will then 
be; 



lA:.s(max) 



= 0.67 X [(7comp - Ksimin)) ' {J-Ks)q\ (11) 



where (J-Ks)o is the intrinsic J-Ks colour of a K or M 

I iant. If we suppose that for a K2 Giant (J-Ks)o - 0.75 
Wainscoat et alJ ll992'). and knowing that in the 2MASS ob- 
servations 7comp varies between nij ~ 12.0 and 15.8 we can 
calculate that the AKs(m-dx) varies between ~ 1.4 and ~ 3.75 . 
From equation Eq. ^\ it is obvious that the highest extinction 
detectable will be limited by the completeness in the / band 
alone. 



For a colour difference between 2MASS and the Galactic 
model to be significant, the difference must be greater than 
the uncertainty in the (J-Kg) colour index in the observations. 
This sets a limit on the minimum extinction detectable. The 
photometric errors in the 2MASS observations depend on the 
Galactic region observed as well as on the local atmospheric 
conditions on the night the observations were made. An uncer- 
tainty of 0.05 magnitudes for the (J-Ks) colour index is a suit- 
able guide, which coiTesponds to a minimum K, band extinc- 
tion of ~0.03 magnitudes (Eq.|3. This error becomes relatively 
small for areas of high extinction, which makes our method 
ideal for extinction studies in the Galactic Plane. 

6.4.2. Maximum and minimum distances 

In order to allocate a distance to a bin of stars from the 2MASS 
survey, we make use of the colour distance relation discussed 
in ^4.11 To do this we need to exclude the dwarf stars from our 
pipeline. As these stars are local, we unfortunately lose most of 
the information on the first kiloparsec. 

The maximum distance varies according to the J and K, 
completeness and column density along the line of sight. 
Generally, we are able to detect the extinction to distances 
greater than 10 kpc, and for many lines of sight well beyond 
this, but for lines of sight with high column densities and low 
completeness (central regions, spiral arm tangents) we may 
only obtain information out to ~ 7 kpc. We also notice that 
lines of sight which contain small number of stars (large abso- 
lute longitudes) contain extinction information to rather short 
distances. Deeper counts with more stars would be needed to 
complete those fines of sight. 

7. Conclusion 

We present the first results of our extinction model in the 
Galactic plane. The Galactic model provides an excellent tool 
for extracting the reddening information locked in the 2MASS 
data. 

Our method does not give much information on the extinc- 
tion in the first kiloparsec from the Sun but does give the large 
scale distribution of the extinction in our Galaxy. Furthermore, 
the J-Ks index is not sensitive to very low extinction zones, 
such as high latitude lines of sight. However, the method can 
easily be modified to use other colour indices, which would 
allow us to adapt the sensitivity to a given region or set of ob- 
servations. An additional restriction is the number of stars in 
the observations for our chosen pixel size. The map could thus 
be expanded to higher Galactic latitudes by lowering the spatial 
resolution. 

Many structures have been identified in the resulting maps, 
including: several spiral arm tangents, dust in the Galactic bar, 
the molecular ring as well as local features such as the Aquila 
Rift. We have given quantitative results for three of these struc- 
tures : 

1 . The dust distribution is found to be asymmetrically warped, 
in agreement with CO and HI observations of the ISM but 
not as pronounced: for positive longitudes the angle is 89° , 
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it starts at 8.7 kpc from the Galactic center and grows with 
a slope of 0.14 , while at negative longitudes, the angle of 
the maximum is at 272° , the starting radius 7.8 kpc and 
the slope 0.11. Hence we confirm that the warp is seen in 
all Galactic components, although the amplitudes appear to 
vary from one component to the other. 

2. By converting our extinction map between 4 and 8.5 kpc 
to a map of hydrogen density, we are able to calculate the 
mean scale height for the interstellar matter in the inner 
Galaxy. We find the scale height of the disc to be 125^'^ 
pc in the region 4 < R < 8.5 kpc. This value of the scale 
height is in agreement with lPrimmel et al.l 112003!) who find 
a value of 134.4 + 8.5 pc. However the dust scale height 
appears to be slightl y smaller than the H I layer , as seen by 
lMa]hotrala995i) and lNakanishi & Sofud ( Eool who find a 
mean value of ~ 160 pc for the same region. 

3. The bulge region is found to contain little absorbing matter 
apart from an elongated structure 5.2 kpc long with an angle 
of 30 + 5° relative to the Sun-Galactic centre direction. This 
dusty bar has a greater inclination than many measures of 
the angle for the stellar bar which may indicate the presence 
of dust lanes preceding the stellar bar at positive longitudes 
and traihng it at negative longitudes. 

In the near future, this 3D map of the location of obscur- 
ing dust in the Milky Way will enable us to constrain dust pa- 
rameters (temperature, size) by comparing the predicted far in- 
frared emission with observations. In addition, this map of the 
Galactic ISM enables us to further constrain the Galactic struc- 
ture parameters via the Galactic model. In particular, the spiral 
structure, the Galactic flare and warp, as well as the size and 
orientation of the bulge and stellar bar will be easier to model 
with this improved estimation of the three dimensional distri- 
bution of interstellar extinction. 
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